Dynamical phenomena in Fibonacci semiconductor superlattices 
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We present a detailed study of the dynamics of electronic wavepackets in Fibonacci semiconductor 
superlattices, both in flat band conditions and subject to homogeneous electric fields perpendicular 
to the layers. Coherent propagation of electrons is described by means of a scalar Hamiltonian 
using the effective-mass approximation. We have found that an initial Gaussian wavepacket is 
filtered selectively when passing through the superlattice. This means that only those components 
of the wavepacket whose wavenumber belong to allowed subminibands of the fractal-like energy 
spectrum can propagate over the entire superlattice. The Fourier pattern of the transmitted part 
of the wavepacket presents clear evidences of fractality reproducing those of the underlying energy 
spectrum. This phenomenon persists even in the presence of unintentional disorder due to growth 
imperfections. Finally, we have demonstrated that periodic coherent-field induced oscillations (Bloch 
oscillations), which we are able to observe in our simulations of periodic superlattices, are replaced 
in Fibonacci superlattices by more complex oscillations displaying quasiperiodic signatures, thus 
sheding more light onto the very peculiar nature of the electronic states in these systems. 

PACS number(s): 72.10.-d, 72.15.Rn, 73.20.Dx 



I. INTRODUCTION 

Since the fabrication of aperiodic semiconductor sru 
perlattices (SLs)_ | arranged according to the FibonacciEl 
and Thue-Morseo sequences, there has been a grow- 
ing interest jru their electronic n grp»erties, both from 
experimentaflia and theoreticalQuO viewpoints. One 
of the most appealing motivation for these studies is 
the theoretical prediction that ideal aperiodic SLs should 
exhibit a highly-fragmented-|electronic spectrum display- 
ing self-similar patterns .tall 3 ] And, in fact, photolumines- 
cence excitation spectroscopy at low temperature reveals 
the existence of a fragmented density of states consistent 
with theoretical predictions 

Other motivation for the study of dynamical phenom- 
ena in aperiodic systems is the following. Electron states 
in periodic SLs spread uniformly over the whole SL 
(Bloch states) and the energy spectrum is composed by 
minibands and minigaps. In the absence of electric field, 
these extended states are characterized by a transmission 
probability very close to unity. When a homogeneous 
electric field is applied perpendicular to the layer plane, 
electronic states become localized (Stark- Wannier states) 
and the energy spectrum consists of equally spaced levels 
(Stark- Wannier ladder). From the perspective of quan- 
tum evolution, Bloch states leadr-to-Bloch oscillations 
(BOs) when the electric is applied. ESEJ The BO time pe- 



riod of the .electronic motion in real as well as in k space 
is given bylB 



where d is the SL constant. High-quality SLs make it 
possible to obtain BO time periods larger than the scat- 
tering time for reasonable values of the applied electric 
field F. Reports of unambiguous experimental evidence 
for BOs in periodic-CaAs-Gai ^Al^As SLs have been re- 
cently appeared, EjO using an experimental method pre- 
viously proposed by von Plessen and Thomas J13 This pic- 
ture is assumed to be no longer valid in Fibonacci SLs 
(FSLs) because in the thermodynamical limit electron 
states are critical instead of extended in flat band con- 
ditions, from a strict mathematical point of view. How- 
ever, since actual FSLs are of finite size, one could ex- 
pect similar transport properties (i. e., high transmission 
coefficient) to that shown by extended electronic states 
since they spread over the whole FSL, although we insist 
that they are not Bloch states. Hence the question as to 
whether BOs will be observed or not in FSL arises quite 
naturally. 

Therefore the aim of this work is twofold. In the first 
place, we provide a complete characterization of elec- 
tronic states in FSLs, giving a detailed description of 
dynamical phenomena of electronic wavepackets, which, 
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as far as we know, have not been reported in the litera- 
ture. In this way, we are led to the conclusion that FSLs 
act as efficient electronic filters, an appealing property 
in order to use them in actual devices of technological 
interest. In the second place, we investigate the possibil- 
ity to observe BOs in FSLs, suggesting the convenience 
of generalizing the concept of periodic BOs to the case 
of quasiperiodic oscillations in order to properly describe 
the dynamical behavior of critical states under the action 
of homogeneous electric fields. 

To this end, in this paper we address the study of time- 
dependent effects in FSLs by solving numerically the 
effective-mass equation for the envelope-function, both 
in the absence of external fields and under homogeneous 
electric fields. To be specific, we consider the problem 
of quantum evolution of electronic wavepackets initially 
localized in space impinging on the FSL. The wavepacket 
dynamics will be properly described by means of the 
time-dependent transmission probability. The transmit- 
ted portion of the wavepacket will be characterized by its 
Fourier transform, aiming to search for particular signa- 
tures arising from the scattering event. To get an estima- 
tion of the spreading of the wavepacket as a function of 
time, we will use the time-dependent inverse participa- 
tion ratio (IPR) as well as the mean-square displacement. 
Finally, since unintentional imperfections appear during 
growth in actual FSLs, we have analyzed a modified ver- 
sion of our model to investigate the possible existence of 
a competition between the long-range quasiperiodic or- 
der and the short-range disorder which could be detected 
by our time analysis. 



II. MODEL 

We consider quantum-well based GaAs-Gai_j,Al a As 
SLs with the same barrier thickness b in the whole sam- 
ple. The thickness of each quantum-well is Ax n — b = 
x n — x n -\ — b, x n being the position of the center of the 
nth barrier and x the growth direction. We will focus 
on electronic states close to the bandgap with = 
and neglect nonparabolicity effects hereafter, so that the 
Ben Daniel-Duke Hamiltonian suffices to describe those 
states. The envelope-functions for electron wavepackets 
satisfy the following quantum-evolution equation 



The particular kind of SLs subject of this work, FSLs, 
can be grown starting from two basic buildingj-blocks A 
and A' by means of molecular beam epitaxy.EJ Here A 
{A') consists of a quantum-well of thickness a (a!) and a 
barrier of thickness 6. The Fibonacci sequence S n is gen- 
erated by appending the n — 2 sequence to the n — 1 one, 
i.e., S n — S n S n -i- This construction algorithm requires 
initial conditions which are chosen to be So = A and 
Si = A' . In this way, finite and self-similar quasiperiodic 
SLs are obtained by n successive applications of these 
rules leading to the ordering A A' AAA 1 A. . . , contain- 
ing N = F n barriers The Fibonacci numbers are gener- 
ated from the recurrence law F n = F n - \ +F n -2, starting 
with Fq — F\ = 1. A few blocks of the resulting SL 
potential Vsl are shown in Fig. |l|. 

Unintentional disorder appearing during growth in ac- 
tual SLs depends critically on the growth conditions and 
it is unknown in most cases. Therefore, one is forced to 
develop a simple model, making reasonable assumptions 
on the type of disorder for each particular sample. For 
instance, islands protruding from one semiconductor into 
the other cause in-plane disorder and break translational 
invariance parallel to the heterojunction. If the in-plane 
average size of these protrusions is much larger than the 
mean-free-path, then capiers only see an ensemble of dif- 
ferent layer thicknesses .113 We model local excess or de- 
fect of monolayers by allowing Ax n to fluctuate uniformly 
around the nominal values a + b or a' + b. For dcfinitcness 
we take Ax n — a(l + We n ) + b or Ax n = a'(l + We n ) + b, 
where W is a positive parameter measuring the maxi- 
mum fluctuation and e n 's are distributed according to a 
uniform probability distribution P(e n ) = 1 if |e n | < 1/2 
and zero otherwise. Note that e n is a random variable, 
even when the mean values of Ax n follow the Fibonacci 
sequence. 



III. NUMERICAL ANALYSIS 

We study the quantum dynamics of an initial Gaussian 
wavepacket 
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The time-independent Hamiltonian Tl(x) is given by 
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where Vsl(^) is the SL potential under flat bands condi- 
tion and F is the electric field. The height of the barrier 
for electrons is given by the conduction-band offset at the 
interfaces AE C . We take the origin of electron energies 
at the GaAs conduction-band edge. 



impinging on the FSL, where the mean kinetic energy 
is (E) = 7l 2 fcg/2m* and AX measures the width of the 
electron wavepacket. The solution of Eq. (0) is given by 



y{x,t) = exp ( -H{x)t ) *(x,0). 



(5) 



The finite difference representation of the exponentia 
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where St is the time step, brings a powerful and high- 
accurate numerical method. In addition, it ensures prob- 
ability conservation,^ which has been used at every time 
step as a first test of the accuracy of results. Boundary 
conditions read $(oo,t) = oo,i) = and we have 
chosen the FSL length sufficiently large to be sure than 
the wavepacket never comes close to the boundaries. 

Transmission properties of the electronic wavepacket 
can be successfully analyzed by means of the time- 
dependent transmission probability Pr(t), which is noth- 
ing but the probability that at time t the electron is found 
to have crossed the whole SL, 

/>oo 

P T (t)= dx\*(x,t)\ 2 , (7) 

where L is the length of the system. In addition, to 
get a complete characterization of the motion of the 
wavepacket, we use the time-dependent inverse partici- 
pation ratio, IPR(i), and the mean-square displacement, 
a{t). The IPR is defined as 

/CO 
dx\V(x,t)\ 4 , (8) 
-oo 

and it gives an estimation of the spatial extent and the 
degree of localization of electronic wavepackets, which 
can indeed provide very much information. Delocalized 
states are expected to present small IPR (in the ballis- 
tic limit, without applied field, it vanishes as while 
localized states have larger IPR. The mean-square dis- 
placement describes how quantum diffusion of wavepack- 
ets initially located in the middle of the FSL takes place. 
The mean-square displacement a(t) is defined as 

/oo 
{x-xf \ip(x,t)\ 2 dx. (9) 
-oo 

with 

/oo 
x \4>(x, t)\ 2 dx. 
-oo 

In the asymptotic regime (t — * oo) one expects <J 2 (t) ~ 
i 7 . The exponent is < 7 < 1 for localized states, 7=1 
for ordinary diffusion, 1 < 7 < 2 for super-diffusion, and 
7 = 2 for ballistic regime. 

IV. RESULTS AND DISCUSSIONS 

A. Zero field behavior 

As a typical SL we have chosen a GaAs-Gao.65Alo.35As 
structure, for which the conduction-band offset is AE C = 
0.25 eV and the effective-mass is m* = 0.067m, m be- 
ing the free electron mass. In our computations we have 
taken a = b = 32 A and a' = 26 A. With our chosen pa- 
rameters there exists only one miniband below the barrier 



in periodic SLs, ranging from 0.102eV up to 0.177 eV. 
Thus the miniband width is much larger than the ex- 
citon binding energy, which amounts ~ 0.01 eV in the 
present SLs. This is a relevant fact, for it has been shown 
that electronic localization (Stark- Wannier states) is sup- 
pressed for low miniband widths. Finally, to compare 
with actual SL's, we have considered the parameter W, 
which governs the imperfection magnitude, ranging from 
up to a maximum of 0.05. This value amounts to having 
protrusions thicknesses of half a monolayer on average. 

Figure || collects the results of a typical simulation of 
a wavepacket for a FSL. First, in Fig. ||(a) we show the 
transmission coefficient, r, as a function of energy for a 
perfect (W = 0) FSL with N = 144 welk, computed 
by means of the transfer-matrix formalismE3 The overall 
structure of the energy spectrum is characterized by the 
presence of four main subminibands (notice that trans- 
mission peaks are clustered around energies ~ 0.120, 
0.155, 0.175, and 0.200 eV). An enlarged view of each 
cluster of peaks shows that the fragmentation pattern 
follows a trifurcation scheme in which each cluster splits 
from one to three subclusters upon going to higher-order 
generations of FSL's. This splittuig-.scheme agrees with 
that previously reported in FLSs.BEa The main question 
is to know how this highly-fragmented energy spectrum 
affects the quantum evolution of a wavepacket incident 
on the FSL. 

Figure ^(a) also shows the Fourier transform of an ini- 
tial Gaussian wavepacket at t = with its average kinetic 
energy (E) = 0.160 eV lying in the centermost submini- 
gap and spatial width AX = 200 A. After the trans- 
mitted electron is found to have crossed the whole FSL 
[~ 6ps, see Fig. |](a)], the Fourier transform of the elec- 
tronic wavepacket changes dramatically in perfect FSLs 
(W = 0). Instead of a smooth function, the Fourier 
transform presents a series of marked peaks. Conspicu- 
ously, the energy of these peaks coincides with the higher 
values of the transmission coefficient, thus indicating that 
the FSL acts as an efficient electronic filter. Notice that 
the Fourier transform also displays the same splitting 
pattern than the energy spectrum, being observable even 
the third level of hierarchy in the upper subminiband. 
A physical understanding of this behavior is achieved 
if one considers that the initially localized wavepacket 
can be regarded as a combination of plane waves in a 
continuous band. Since the dispersion relation (energy 
versus wavenumber) is self-similar with a hierarchy of 
splitted subminibands separated by well-defined mini- 
gaps [see the transmission coefficient shown in Fig. ||(a)], 
only those components whose wavenumber belongs to an 
allowed subminiband can propagate over large distances 
and, consequently, contributing to the trasmitted part of 
the wavepacket. What it is most important for practi- 
cal purposes, we have found that unintentional disorder 
does not severely affect filtering properties, as shown in 
Fig. ^(b) for W = 0.05. Although an overall reduction of 
the transmitted components is found, signatures of the 
above mentioned level splitting are still clearly observed 
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in the Fourier patter, particularly at the central energy 
region around 0.155 eV. 

Quantum evolution of electronic wavepackets will de- 
pend upon the system length since the fragmentation of 
the energy spectrum is higher on increasing N. Fig- 
ure [|(a) shows the time-dependent transmission prob- 
ability as a function of time for perfect (W = 0) FSLs of 
various lengths. For comparison, it should be borne in 
mind that the transmission probability vanishes in inten- 
tionally disordered SLs for moderately large sizes, thus 
providing further evidences of the differences between 
random and aperiodic systems. The occurrence of the 
plateau for larger times indicates that the transmitted 
wavepacket has crossed the whole FSL. Thus, the larger 
the FSL length, the later the plateau appears, as was to 
be expected. Interestingly, the value of the asymptotic 
probability decreases on increasing the FLS length. This 
reduction of the transmission probability means that fil- 
tering effects are stronger as the fragmentation of the 
energy spectrum is higher. That is to say, the equivalent 
miniband- width, defined as the sum of all the allowed 
subminibands, decreases as a power of N as ap-eonse- 
quence of the quasiperiodic topology of the FSLO and, 
therefore, more and more components of the wavepacket 
are back reflected. Figure ||(b) presents the results for the 
time-dependent transmission probability as a function of 
time for imperfect (W = 0.05) FSLs of various lengths. 
A comparison with Fig. ||(a) indicates a decrease of the 
transmission probability since short-range disorder tends 
to localize particles, according to the Anderson theory. 
However, this effect is almost unnoticeable for relatively 
short (say N smaller than 89) FSL, in agreement with our 
previous estimations based on the study of the equivalent 
miniband- width. c3 Therefore, we are led to the conclu- 
sion that moderately large fluctuations cannot destroy 
filtering effects of actual FSLs and that they can safely 
be used for practical purposes. The study of the mean- 
square displacement in perfect and imperfect FSLs pro- 
vides additional support in this sense. We have placed 
the initial Gaussian wavepacket in the middle of the FSL, 
with AX = 20 A and (E) = 0.160 eV. In all cases we 
have observed the super-diffusive regime a 2 (t) ~ t 1 with 
7 ~ 1.2, as illustrated in Fig. |j. This result shows that 
the time-dependent properties are quite similar for both 
perfect and imperfect FSLs as large as samples contain- 
ing N = 377 wells. It is also interesting to note that the 
exponent 7 obtained in our study, which is based on a 
realistic continuous model, is lower than that obtained 
by considering tight- binding Hamiltonians£9 



B. Homogeneous field effects 

We now comment on our results when a homogeneous 
electric field is applied perpendicular to the layers. One 
expects that BOs are to be observed in periodic SLs in 
the limit of high electric fields. The localization length of 



Stark-Wannier states, the static counterpart of the BOs, 
is of the order of AE/eF, AE being the width of the al- 
lowed miniband, so that we can roughly estimate that the 
condition AE/eF <~ a + b establish the high field regime. 
In our SL this field amounts F ~ 100 kV/cm. Figure @(a) 
presents the results for the IPR when the initial Gaus- 
sian wavepacket with AX = 20 A and (E) = 0.160 eV 
is located in the middle of the periodic SL, the electric 
field being F = lOOkV/cm. The IPR displays periodic 
oscillations with marked peaks at times tk — fcre, where 
k is any arbitrary, nonnegative integer and tb = 0.065 ps. 
Thus the IPR is bounded below, indicating that the 
wavepacket is localized (dynamical localization) but its 
spatial extent varies periodically in time. Notice that the 
value of the oscillation period is in excellent agreement 
with the theoretical prediction tb = 2nh/(eFd) given in 
(|l|). It is also worth to mention that the initial state 
\l/(x,0) is not completely restored after each oscillation, 
as it should be if interband tunneling were negligible. Efl 
Thus we are led to the conclusion that intraband tunnel- 
ing plays a role in our periodic SL. 

Results corresponding to a FSL with the same initial 
condition as before are also shown in Fig. ||(a). First 
of all, we observe that at short times we can detect two 
oscillations coinciding with the positions of the first two 
BO peaks, but at larger times periodic BOs are com- 
pletely absent in FSLs. The absence of Bloch oscillations 
in FSLs simply reflects the fact that their extended states 
are no longer Bloch states. Bloch states are character- 
ized by a periodic pattern, but this is not the case in 
the FSL, where critical states spreading over the whole 
system show quasiperiodic patterns.EJ In this sense, sig- 
natures of of a quasiperiodic oscillation pattern in the 
IPR corresponding to the FSL can be seen from Fig. ||(a) 
in the interval 0.5 < t < lps. To this end, we intro- 
duce the ratio p/q, where q is the number of oscilla- 
tions of the FSL IPR comprised in a given number p 
of periodic BOs. In this way we obtain the following se- 
quence p/q = 2/3,3/5,2/3,5/8,2/3,7/11,2/3... which 
converges to 2/3, an approximant of the inverse golden 
mean tq = (\/5— l)/2. 

This facts can be described by means of the following 
scenario, based on our previous discussion on the local- 
ization degree of electronic states under an applied elec- 
tric field. For such high electric fields the spatial extent 
of the electronic states is of the order of the SL constant 
and, consequently, the electron cannot see the long-range 
quasiperiodic potential. In fact, differences between pe- 
riodic and Fibonacci SLs at short times are actually very 
small [see Fig. |(a)] . As soon as some components es- 
capes from the localization region (for instance, those 
states above the barrier has a high transmission coeffi- 
cient) the effect of quasiperiodicity takes place. To get 
further confirmation of this assertion we have also stud- 
ied the case of lower electric fields, as shown in Fig. ||(b) 
for F — lOkV/cm, when the localization length is about 
ten SL periods. It becomes apparent that no signatures 
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of BOs can be detected even at short times. Moreover, 
Fig. ||(b) also displays the result for the periodic SL and 
the same electric field. Notice that the initial states is 
completely restored after each BO, the suggesting that 
in this regime interband tunneling is actually negligible. 



lished between quasiperiodic nanoelectronic devices and 
information theorytil pave the way to a very exciting and 
promising line of cross-disciplinary research. 
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V. CONCLUSION 

In this paper we have studied quantum dynamics of 
wavepackets in Fibonacci SLs. After solving numerically 
the time-dependent effective-mass equation arising from 
the Ben Daniel-Duke Hamiltonian, we have found that an 
initial Gaussian wavepacket undergoes complicate scat- 
tering events at the FSL. In particular, the Fourier spec- 
trum of the transmitted part of the wavepacket is no 
longer smooth but presents many marked peaks. These 
peaks correspond to the allowed subminibands of the en- 
ergy spectrum, thus being a clear indication that the FSL 
acts as a fine electronic filter. Transmission properties 
under zero bias have also been successfully characterized 
by means of the transmission probability given by (Q). 
From these results we have learned that unintentional 
disorder arising during growth have no noticeable effects 
on the filtering capabilities of short FSLs, because they 
only cause an overall decrease of the transmitted ampli- 
tude. On the other hand, the spatial degree of localiza- 
tion of wavepackets driven by an electric field has been 
properly described by means of the time-dependent IPR: 
As a check, in periodic SLs we have obtained the dy- 
namical localization fields as well as BO's. On the con- 
trary, quantum dynamics in FSLs also exhibits dynam- 
ical localization although it turns out to be much more 
intricate. In particular, no evidence of periodic Bloch 
oscillations was observed. Instead, the very concept of 
BO should be revisited in order to include the possibility 
of quasiperiodic pattern oscillations emerging from the 
quasiperiodic nature of the system. This result puts an- 
other piece of evidence about the very .peculiar nature of 
electronic states in Fibonacci systems.Ell 

As a concluding remark, we want to stress that with 
this work we have sufficiently demonstrated the existence 
of distinct physical observable consequences of singular 
energy spectra like that of the FSL. For instance, the 
fragmentation of the SL minibands could serve as the 
basis for electronic filters, in a similar fashion to metallic 
Fibonacci multilayers acting as selective filters of soft X- 
ray radiation.EZl It is quite clear that the output of such 
a device is but an image of its fractal spectrum, which, 
in turn, is intimately connected with the quasiperiodic 
nature of the SL and hence with the information it con- 
tains. Therefore, aside from the possibility of building 
the specific filter-like devices we already mentioned, de- 
signed with this or other quasiperiodic sequence accord- 
ing to the desired application, it might as well be that this 
kind of systems can be used in transmitting or process- 
ing information. This and other relations already estab- 
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FIG. 2. (a) Transmission coefficient versus incoming en- 
ergy for a GaAs-Gao.65Alo.35As FSL with N = 144 wells 
(thin solid line). Fourier transforms of the initial wavepacket 
(AX = 200 A and (E) = 0.160 eV, dashed line) and the trans- 
mitted wavepacket for perfect FSL (W = 0, thick solid line) at 
t = 6ps are also shown in arbitrary units, (b) Fourier trans- 
forms of the trasmitted wavepackets for a perfect (W = 0, 
solid line) and imperfect (W = 0.05, dashed line) FSLs at 
t = 6ps. Notice in (a) the perfect coincidence between the 
filtered wavepacket and the allowed minibands. Parameters 
are a — b = 32 A and a' — 26 A. 
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FIG. 1. Schematic diagram of the conduction-band edge of 
the GaAs-Gai_ B AlyAs Fibonacci superlattice. 
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FIG. 3. Transmission probability versus time for (a) per- 
fect W = and (b) imperfect W = 0.05 GaAs-Gao.65Alo.35As 
FSL with various number of wells: From top to bottom 
iV = 34, 55, 89, 144, 233. Other parameters are the same as 
in Fig. 2 
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FIG. 4. Mean-square displacement versus time for perfect 
(W = 0, solid line) and imperfect (W = 0.05, dashed line) 
GaAs-Gao.65Alo.35As FSL with N = 377 wells. Other param- 
eters are the same as in Fig. 2. The mean-square displacement 
grows in time as a power law ~ t 1 . 
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FIG. 5. Inverse participation ratio as function of time 
for an initial Gaussian wavepacket placed in periodic (dashed 
line) and Fibonacci (solid line) GaAs-Gao.65Alo.35 As SLs, 
subject to an electric field (a) F — lOOkV/cm and (b) 
lOkV/cm. Other parameters are the same as in Fig. 2. The 
occurrence of Bloch oscillations in the periodic SL is apparent. 
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